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1,0 INTRODUCTION 


The purpose of this document is to provide the user with a manual 
containing all information necessary to operate the Steady-State Attitude 
Control Propulsion System Computer Program (SSACPS). The SSACPS was 
developed as a general computer simulation model capable of modeling a 
large variety of attitude control propulsion systems (ACPS) • This model 
has been developed under MSC/TRW Task 705-1 along with the ACPS Transient 
Computer Program. 

The SSACPS was developed with the capability to simulate operation of 
ACPS components connected in a large variety of configurations. The coirh- 
ponents are connected within the program by the user, using data required to 
define initial and operating conditions. The component models consist of 
gas generators, combustors (monopropellant and bi-propellant) , lines 
(adiabatic, isothermal, and heat transfer with f riction- liquid /gas) , heat 
exchangers, junctions, regulators, pressure boundaries, and liquid and gas 
valves. The system model uses hydrogen and oxygen properties from detailed 
NBS property subroutines. The data are input through the numerical input 
subroutine TRWLOD and the system equations are solved by the nonlinear 
equation solver, TRW BEGS. The constraint equations for each component set 
and the state variables are selected by the user. Tliis allows a great deal 
of flexibility in solving for either fluid state points or component dim- 
ensions. 

This document is comprised of two volumes and is to be used as a supple- 
ment to the Engineer's Manual, (Reference 1). Volume I contains the user’s 
and programmer’s information and Volume II (Reference 2) contains a copy of the 
program listing. Reference 1 contains more detailed information on each compo- ^ 
nent model. The SSACPS was developed for use on MSC’s UNIVAC 1108 computer. 
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2.0 PEDGRAM DESCRIPTION 

2. 1 Program Definxtion 

The SSACPS is a flexible steady-state attitude control propulsion 
simulation computer program. The program consists of engineering equations 
simulating the fluid states using properties of hydrogen and oxygen. The 
combustor gases are modeled using curve fit data. ACPS components are 
available for combxnation in any type of ACPS simulation. 

The SSACPS program is stored on a PCE tape which contains the 
entire program. The program includes subroutines which model individual 
components, thermodynamic property data, a nonlinear equation solver, and 
a set of driver subroutines. 

TRt^LOD (card input) input is used to enter the component data and 
system configuration. TRl-ILOD is described in Section 4.8 and the descrip- 
tion of the card inputs is presented in Section 3.1.3. 

2. 2 Method of Solution 

The theory and mathematical derivations for this program are des- 
cribed in Reference 1. 
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3. 0 USERS INFORM^^TIOK 

The following sections present detailed descriptions of the types of 
input data and rts format. 

3. 1 Input Description 

3.1.1 Solution of Nonlinear Equations: Input 

Eor normal usage IFIN=0. IFIN^^O means that the finishing procedure 
only is to be executed. The purpose of this procedure is to examine the 
function E in a neighborhood of a given point where, presumably, F(X^) 
is a minimum. This is accomplished by perturbing each component of by 
plus and mxnus small amounts, evaluating the corresponding functions, and 
printing all information as output via the system output tape. The routine 
was designed such that BEGS may be called with IFIK=0, and then immediately 
called again with IFINt^O without changing any other parameters. This 
procedure allows the user to test the solution just obtained. 

The parameters are defined in Section 3.1.1 for the normal usage option 
(1FIN=0) , and in Section 3.1.2 for the finishing procedure. For either 
option the parameters H, BOUNDS, JVECTR, and EPSLIN may be zero upon entry; 
"reasonable" values will then be substituted by the program. 

Printout of information concerning any difficulties encoiantered by the 

program in solving the equations can be obtained by setting EPSLIN<0. 

3. 1.1.1 Solution Normal Usage (IFIH=0) 

STATE = State variables Xx (Xi...Xjjj) used in non linear 

equation solver; first value input is the name of the 
variable, second and third values are the location 
of the variables, and the fourth value is the users 
starting guess. (State output is the program solu- 
tion) 
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■Example: 

STATE = WD0TG, 1,1,5 . ,AKEAT,1,1, ,3 
flowrate WDJfTG, line 1, node 1, 1st guess 
5.1bm/sec, area throat, combustor 1, Dunui^ 
value 1, 1st guess. 3 square inches. 


The following 

list presents the state variables used in the 

non-linear 

equation solver: 



WDOTG(II,JJ) 

Elowrate, Line II, node JJ 

Ibm/ sec 

MACH (II, JJ) 

Mach nxmber, line II, node JJ 

- 

PG(II,JJ) 

Pressure, line II, node JJ 

Ibm/in^ 

XK(I) 

Area of value I times discharge coefficient 

in^ 

PRE(I) 

Discharge pressure of regulation I 

2 

Ibm/ in 

DIALI(II) 

Diameter , line II 

CM 

TTEMP(II,JJ) 

Temperature, line II, node JJ 

deg-R 

TTTMP(II,JJ) 

Total temperature, line II, node JJ 

deg-R 

PGPBL(I) 

Pressure at Pressure Boundary I 

Ibm/ in^ 

ARINJ(I,II) 

Area Injector, Combustor I, line II 

. 2 
in 

PCI (I) 

Combustor (I), Chamber pressure 

Ibm/ in^ 

AREAT(I) 

Area throat Combustor I 

, 2 
■in 

THRT(I) 

Thrust Combustor I 

Ibf 

PDELP (I) 

Pressure rise through turbopump (I) 

psia 

LENGTH (I) 

Length heat exchange tubes 

in 

RI(I) 

Inner radius of concentric heat exchanger 
tubes 

in 

R0(I) 

Outer radius of concentric heat exchanger 
tubes 

in 

HR(I) 

Hydraulic radius 

in 

P(I) 

Accumulator Pressure 

Ibm/in^ 


(h2j...,hj{) [input] is a vector of perturbations on X 
for the initial functxon evaluations. If h^ = 0, then 
hi = ,01(xi')will be automatically substituted; if Xi = 0, 
hi = .001 will be substituted Instead. 

H [output] is set-to (0,0,.,., 0). 
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IBNDS 


B0BNDS 

B0UNDS 

M 

N 

MAXIT 

IDER 

JVECTR 


EPSLN 


BL0CK 


ITERS 


[input] is zero if the solution to the linear least squares 
solution by (LEGS) is to be unconstrained. 

IBNDS is nonzero if the least square solution is to be con- 
strained by the use of B0UNDS (described below) . The recom- 
mended usage here is IBNDS^O, for the routine will evaluate 
B0UNDS and if too small, will increase its values to a size 
where they have no effect. 

(bi,...»by) [input] IS the vector of constraints for 
bounding the least squares solution (IBNDS). If any 
bj^=0, then bj^=10 will be automatically substituted. 

[output] is set to (0,0,..., 0). 

[input] xs the ntimber of variables. 

[input] is the number of functions. 

[input] xs the maximum number of iterations allowable. 

[input] is zero for the secant method. IDER is non- 
zero if the user computes the matrix of partial 
derivatxves in Component Sub ro untines. i 

[input] is a vector of integers, which allow 
the user to hold fixed a particular combination of the 
variables. If x^ will be disregarded in the 

iterative procedure, but will be present as input to 
FC0MP. If ji=0, x£ is taken to be a true variable and 
is solved for. 

[input] is a convergence test. The routine will exit 
if the following test xs passed: 

[i Xi I < e I Xj_ I , i “ 1,M 

, f I epSLnI if EPSLN 4 0 

where s = \ S 

( 10“ if EPSLN = 0 

Here A xi is the newly computed correction to the current 
i-th component of the xterant. If EPSLN < 0 a message is 
printed out to indicate an unsuccessful secant iteration. 
See Step E in the Appendix. 

is a FORTRAN vector used for temporary storage. The 
DIMENSION of this vector must be at least (3M +4) (M + 1) 

+ -N (2M + 1) . 

[output] the magnitude of ITERS is the actual number of 
iterations completed before returning to the calling 
program. If ITERS = MAXIT, the "solution" should be re- 
garded with suspiexon. If ITERS is negative, iteration was 
terminated due to inability to find suitable function 
values . 
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2 ^ 

CONSTRAINTS/F (fl, . . . ,£n) [output] is ?(X) such that jlF]] is a 
minimum . 

FN^BM [output] is the sum of the squares of 

NOTE: If the process didn't converge, X and F correspond 

to the smallest value that was obtained for FN^EM. 


3. 1.1. 2 Solution Using Finishing Technique (IFIN 4 0) 


STATE [output] is a trial solution to be tested by the finish— 

ing procedure. _ If the perturbed state variable produces 
the best solution, the perturbed varxable is used as the 
state variable for the next pass. If no better solution 
was found the original state variable is used. 


H [input] is a vector of small perturbations. If hi = 0, 

then the value x^ will be substituted (see EPSLN, Section A. 

H [output] - (0,0,..., 0). 


M 


See Section 3.1.1 


N 

See Section 3.1.1 

JVECTR 

See Section 3* 1.1 (one may bypass perturbing any of the 
variables) . 

EPSLN 

See H [input2 above. 

F 

[output] is F(X) 

FN0RM 

iB IlFll' 

3.1.2 General 

Input 


The SSACPS input is used in a way which defines the selected system. 
Each component is defined in the SYSCOM array. This array defines the 
order in which the program calculations will proceed. For example, calling 
SYSC0M=TAHK, 1, PIPLL,4, PEL, 3 would set up an ordered three- component 
connection of a tank, liquid line and a constant pressure boundary. Each 
component must have all of the inlet connections defined by input prior 
to calling the component in the array SYSCOM. For example, before calling 
combustor 1, the following system setup must be set by input: SYSTEM = 
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PIPLjl, PXPL,2, CHAM,1* This represents the connection of pipes 1 and 2 
to combustor 1. Prior to calling a turbopump the following information 
must be set by input; SYSC0M=PIPLL,1, PIPL,2, PXPL,3, TUPBOP,!, PIPL^4. 
This turbopump input identifies a liquid line PIPLL at the inlet to the 
pump and PIPL 2 and 3 as the gas lines input to the gas generator (GG, 
modeled as one of the combustors^ CHAM). Note that the GG is automatically 
connected to the turbine inlet of the turbopump. The order of the com- 
ponent selection is also the order which the calculation in the program 
will proceed. Both the individual components and their connections must 
be defined in detail. 

The required inputs for each component are defined in Sections 
3.1.4 through 3.1.15. The state variable (STATE) ^ and functional con- 
straints (PUNCT), defined in Section 3.1*2 are used in the program input. 
The state variables (STATE) are the unknowns in the system and the func- 
tional constraints (FUNCT) are the constraint equations used by the non- 
linear equation solver to obtain a solution. An example of the input 
used to define the state variables are as follows: 

STATE - AEINJ,1,1,-9, ARINJ,1,2,.8 AREAT,1,1.3. 

The first value ARINJ represents the area of the injector in chamber 1, 

2 

line 1 and the first guess used in the program value is 0.9 inches . 

ARINJ is the area of the injector connected to chamber 1 and line 2. 

AEEAT is the area of the throat of chamber 1, the second value of 1 is a 

dimmy value necessary for input only and the initial guess for the throat 
2 

area is 3. inches . The associated functional constraints to go with the 
state variables defined above as follows: 

FUNCT - C0NCT,1, CONCT,2, C0NCT,3 
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The first constraint, CONGT 1, is the injector area constraint for line 1 
in the first combustor; COKCT 2, is the injector area constraint for line 2 
combustor 1; and CONCT 3 is the throat area, chamber pressure or thrust 
constraints for combustor 1, Proper selection requires that the number of 
constraints match the number of state variables, A sample input case is 
presented in Section 3,2 

The following represents the description of data cards for TR\^LOD 
input package.. 

Input of Numerical Data may assume the following foinns in columns 
1-72; 

(a) User defined argument name == constant 

(b) User defined array name = set of constants (separated by 
ccjmmas ) 

(c) User defined argument = subscript, set of constants 
(separated by commas) or 

User defined argument, subscript = set of constants 
(separated by commas) 

(d) User defined arra^^ name = dimension of variable, 
subscript, set of constants (separated by commas). 

The dimension should appear bet^^/een slashes and indicate 
the number of rows in a t\^o- dimensional array or the 
number of rows and columns in a three-dimensional array* 

The dimension and subscripts must follow the equal sign. 

(e) User defined array name nRV (value V to be repeated 
n times) 

(f) User defined argument name = holler! th characters 
(alphanumeric only, first alphabetic) 

(g) User defined argument name = 'hollerith characters’ 

(limited by apostrophes) 

(h) User defined argument name - octal constant 


10 



Any combination of the eight types of data input described above may be 
used in a data set. Data items must be separated by commas* The end of 
a data set is signified by ENDCAS. An exponential input format may be 
used, e.g., A = 1, E-15, If the decimal point is omitted and E, D, or 
sign appears in the value the number will be converted, assuming the 
decimal point is to the right of the last numerical character (in the 
base number). An absence of a decimal point, an E, a D, or an imbedded 
sign signifies an integer value. An asterisk following the data causes 
scanning of that data card to stop, permitting comments to be placed 
after the value. The scan of a card also terminates if an error is 
detected. 

Data may be represented on Numerical Data cards in any of the follow- 
ing forms: 

lA 
X 

X(2) 

D 

D 

DX 
DX 
DX 
A 

integer, double precision, hollerith, or 
octal) 


= 123* THIS IS A COMMENT 
= 1.23E1,2.9 
=2.9, (1) 123-1 

= /5/(2,3), 17.1* (where D is dimensioned (5,n), and 
17.1 is to be stored in D(2,3)) 

= /5 ,3/ (2,2,3) , 17.1* (where D is dimensioned (5,3,n) 
and 17.1 is to be stored in D(2,2,3)) 

= 1.23D4 
= 123D2 
= 1.23D+04 

= 3R14.1, 4R2* (fields to be repreated may be real. 
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ICON = 0040045000000 

ICON = 0040045* (twelve digit word is filled with six zeros 

to the right) 

A = ABCD* (one word is filled with the four characters 

follo^jed by t^^o blanks) 

A =* *bb(ABCD)bb’ (two \fords as follows; bb(ABS and D)bbbb) 


FUNCT 

SYSC0M 


NSYS 


Functions array (2,N) Stored column-wise as name and 
subscript 

Component configuration input. The first value is the 
name of the system component, the second value is the 
number of the component. 

Total number of system components in Syscom. 


3.1.3 System Components. The following is a list of system component 
modules available in this program for compilation into any desired attitude 
control propulsion system; 


PBL 

Pressure Boundairy (6 press boundaries) 


PIFL 

Gas pipe line with friction and heat transfer (30 
lines maximtim including PIPLL, ADIAB, ISOTH lines) 

PIPLL 

Liquid pipe line 


ADIAB 

Gas adiabatic pipe 


ISOTH 

Gas isothermal pipe 


JUNGL 

Junction (6 junctions) 


CHAM 

Bi propellant combustor or gas generator (6 
including (CMON) 

combustors) 

CMON 

Mono propellant combustor 


TANKD 

Accumulator or storage tank (5 tanks) 


TUPBOP 

Turbopump (5 turbopumps) 


VALVL 

Liquid valve (20 valves including VALVG) 


VALVG 

Gas valve 


PEEG 

Pressure regulator (10 regs) 


FEEDB 

Feedback constraint subroutine (5 values) 


HEATX 

Heat Exchanger (5 each) 
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3.1.4 Constraints 


The constraints (F) are incorporated in the system component models 
and are utilized by the nonlinear equation solver to obtain the smallest 
constraint such that | | F | j 2 is a minimum. 


Symbol 

Description 


Dimensions 

PB(I) 

Pressure boundary constraint 


6 

JUN(I) 

Flowrate at junction const raint 


6 

C0NCT(I,J) 

Combustor oxidizer flowrate constraint (IjJ) 
Combustor fuel flowrate constraint (2,J) 
Combustor chamber pressure constraint (3,J) 

3,6 

POWC(I) 

Turbopump power balance constraint 
(pump power = turbine power) 


5 

WTGC(I) 

Turbxne flowrate constraint 
(turbine flow = gas generator flow) 


5 

FBWC(I) 

Feedback flowrate constraint 


6 

FBPC(I) 

Feedback pressure constraint 


6 

FBTC(I) 

Feedback temperature constraint 


6 

TH0C(I) 

Heat exchanger outlet temperature 
constraint 


5 

CGGTC(I) 

Chamber (GG) combustion temeprature 
constraint to heat exchanger 


6 

CGTMK(I) 

Gas inlet pressure constraint to tank 

5 

CPJU(I,J) 

Junction inlet lines pressure constraints 

10,5 

CLTANK(I) 

Inlet liquid pressure constraint to 

tank 

5 

3.1.5 PBL Pressure Boundary (Reference 1, Section 3 

.4.3.2) 


Symbol 

Description 

Units 


PBL(I) 

System component name for pressure 
boundary I 


lEND(I) 

END to which pressure boundary I is 
connected (l-upstream, 2-downstream) 


IPB(I) 

Tells which line is connected to 
pressure boundary I 

6 
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Symbol 

Description 

Units 

Dim. 

PGPBL(I) 

Pressure at pressure boundary 

PSIA 

6 

TT(I) 

Switch for total temperature input 
to PBM (=1 Input total temperature 
=2 Input fluid temperature 


6 

TTTEMP (II, JJ) 

Static temperature 


20,30 

TTTEMP (II, JJ) 

Total temperature 

°R 

20,30 

PB(I) 

Pressure boundary downstream 
constraint 


6 

3.1t6 Junction 

input (Reference 1, Section 3.4.7) 



Symbol 

Description 

Dimensions 

JUNCL(I) 

System component name for junction I 

10 

Jira(i) 

Flowrate constraint in Junction subroutine 10 

NLINJU(I) 

Number of pipes at the Ith junction 

10 

IJUNC(I,J) 

Array which denotes connection to junction I 10,6 

of J lines, upstream lines must be input 
first positive and the downstream lines 
negative 

3.1.7 Valve/Orifice (Reference 1, Section 3.4.5) 



Symbol 

Description 

Units 

Dim. 

VALVG(I) 

System component name gas valve I 


20 

VALVL(l) 

System component name liquid valve I 


20 

VAD(I) 

Line to which domstream end of valve 
is connected 


20 

VAI (I) 

Line to which upstream end of valve 
is connected 


20 

XK(I) 

Discharge coefficient times area of valve 

T 2 

In 

20 

3.1.8 Mono /Bi-Propellant Combustors (Reference 1, Section 

3.4.11) 


Symbol 

Description 

Units 

Dim. 

AEEAT(I) 

Chamber Throat Area 

1N**2 

6 

ARINJ (I) 

Area Injector 

IN** 2 

6,2 

PCI (2) 

Chamber Pressure constraint 

PSIA 

6 
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Units Dimensions 


Symbol 

AJ(I) 

PDEDP(I) 

NCHAM (I) 
ICHAM(I,N) 

CV(I,J) 

THRT 

PATM 

CMON(I) 

CTW 


LAMDA(I) 
GDC (I,J) 
ISPI(I) 
AREAE(I) 

CHAM (I) 
NCOEF 

MW2 (I) 
TC2(I) 


Description 


Switch for input PDROP (AJ=L, 6 

Input PDROP) 

Pressure drop across injector PSIA 6 

PSIA 

Number of chambers 6 

Number of lines connected to the 6,2 

nth inlet of combustor I 

Injector discharge coefficient 6,2 

Engine Thrust constraint Ibf 6 

Atmospheric pressure PSIA 1 


*^Mono-Propellant Only** 
System component name for mono- 
propellant combustor I 
Constraint equation switch for 
combustor' 

=1 Specific xmpulse constraint 
=2 Chamber measure constraint 
=3 Thrust constraint) 


Divergence nozzle coefficient 6 

Chamber flowrate discharge coefficient 6,2 

Specific impulse constraint sec 6 

Area exit chamber in**2 6 


**Bi-Propellant Only** 

System component name for bi-propellant 
combustor I 

Number of coefficients in curve fits 1 

for propulsion parameters 

Molecular weight coefficients in Deg.R 4 

curve fit equation combustor 

Temperature coefficients in curve 4 

fit equation 
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Svmbol 

Description 

Units 

Dim. 

CP2(I) 

Specific heat coefficients in curve 
fit equation 


4 

CS2(I) 

CSTAR coefficients in curve fit equation 

ft/sec 

4 

IS2(I) 

ISP coefficients in curve fit equation 

sec 

4 

C0NCT(1,I) 

Thruster (I) flowrate, w , constraint 
equatxon 

sec 

3,6 

C0NCT(2,I) 

Thruster (I) flowrate, w^, constraint 

sec 

3,6 


equation 



CONCT(3,I) 

Pc, thrust constraint equation 

sec 

3,6 

3.1.9 PIPL Gas 

Lxne (Reference 1, Section 3, 4*1) 



Symbol 

Description 

Units 

Dim. 

ADIAB(I) 

Sys.tem component name adiabatic line I 


30 

PIPL (I) 

System component name gas line I 


30 

PIPLL(I) 

System component name liqiiid line I I 


30 

ISOTH(I) 

System component name isothermal line I 


30 

ALPHA (II) 

Angle line (11) makes with 
horizontal 

Peg. 

30 

TH(II) 

Wall thickness of line (II) 

n* 

4.20, 

KPIPL 

Number of lines 


1 

DIALO (II) 

Outside diameter line II 

n. 

30 

GEE 

Acceleration due to gravity 

2 

ft/ sec 

1 

HNPR 

Exponent on Prandtl number for line 
inner wall film coefficient 


1 

HKRE 

Exponent on Reynolds number for line 
inner wall film coefficient equation 


1 

HOEX 

Exponent on Prandtl+Grashof numbers for 
line outer film coefficient (stored 0.25) 

1 

TENV 

WTFCON(I) 

Temperature of environment 

Wall thermal conductivity ^ 

Deg. R 
BTU 

1 

30 

Et-hr°R 

XLENGL(I) 

Length of the Ith line 

In. 

30 
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Symbol 


Description 


Units Dimensions 


MACH (II, JJ) 

Mach number line II node JJ 

20,30 

NPLINE(II) 

Liquid in line II (1-OX, 2s=Fuel) 

30 

HRS (I) 

Switch for hydraulic resistance for 
Ith line 1=0 program calculates 
HR 

=1 input HR(liquid line) 

30 

HR(I) 

2 2 

Hydraulic resistance (liquid line) Sec /ft-in 

of Ith line 

30 

laxd) 

Total loss coefficients for bends, 
contractions, etc. (liquid line) 

30 

3. 1,10 Heat Exchangers (Reference Section 3* 4 ,15) 


i-> 

HE ATX (I) 

System component name for heat exchanger 

5 

TYPE (I) 

Type of heat exchanger, (3= 
crossflow, 2= parallel flow, 
1= counter flow) 

5 

LI (I), L0(I) 

Inner and outer fluid number 
(1= hydrogen) 

(2= oxygen) 

(3= ^ 

5 

XT, XO 

Inner and outer fluid inlet quality 

5 

HE0I(I) 

Number of the line connected to outer line 
of the Ith* heat exchanger inlet 

5 

HEII (I) 

Number of the line connected to inner line 
of the Ith* heat exchanger Inlet 

5 

HE00(I) 

Number of the Ixne connected to outer line 
'of the Ith* heat exchanger outlet 

5 

HEI0(I) 

Number of the line connected to inner line 
of the Ith* heat exchanger outlet 

5 

HGII (I) 

Number of combustors connected to Ith 
heat exchanger inside line 

5 

BGI0(I) 

Number of combustors connected to Ith 
heat exchanger outside line 

5 

TH0C(I) 

Temperature constraint for heat exchanger 
cold side discharge temperature 

5 


* If the heat exchanger is flowing into the line, the line number 
is input as a positive value, if the flowrate is going from t e 
line into the heat exchanger, the flowrate is input as a 
negative value. 
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**Coimter Flex'? Heat Exchanger** 


Symbol 

Description 

Units 

Dim. 

P(l) 

Pitch distance betx^een tubes 
center to center 

In 

5 

ARFL(I) 

Flow area cross section of heat 
exchanger outer tube 

in2 

5 

NT (I) 

Number of tube wraps 


5 

DTB(I) 

1 

Diameter of tube bundle 

In 

5 

HRS (I) 

Switch for hydraulic input 
radius (=1 input, =0 calculated 
in program) 


5 

HR(I) 

Hydraulic radius 


5 

KLT(I) 

Loss for entrance, bends, and 
discharge coefficient 


5 

IT (I) 

Length of inner tube 

In 

5 

LS(I) 

Heat exchanger shell length 

In 

5 

DHYO(I) 

Hydraulic dia for crossflow 
helicoil heat exchanger 


5 

DSI (I) 

Diameter of heat exchanger shell (inside) 
‘ total heat transfer for crossflow heat 
exchanger 



**Parallel or Crossflow Heat Exchanger** 



Symbol 
LENGTH (I) 

~ Description — 

Heat exchanger length In 

Dimensions 

5 

EL (I) 

Inner tube radius 

In 

5 

RO(I) 

Outer tube radius 

In 

5 

THICK(I) 

Wall thickness of inner tube 

In 

.5 

3.1.11 Tank 

(Reference 1, Section 3.4.3) 



TANKD(I) 

TIMEND 

Svstem coinDonent name for 
accumulators/tanks I 

Time to end program 

sec 

1 
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Symbol 

Description 

Units Dimensions 

BELT 

Time increment 

1.0 sec 

1 

V7D0TG(II,JJ) 

Inlet liquid flowrate 

Ibm/sec 

30,20 

WD0TG(II,JJ) 

LDMI(I) 

Outlet liquid flowrate 

Number of lines, connected to tank 
gas inlet’ 

Ibm/sec 

30,20 

5 

LDMO(I) 

Number of lines, connected to tank 
gas outlet 


5 

LDWI(I) 

Number of lines, connected to tank 
liquid inlet 


5 

LDWO(I) 

Number of lines connected to tank 
liquid outlet 


5 

AVElfr(I) 

Relief valve flow area 

. 2 
in 

5 

CPI (I) 

Specific heat -- species 1 

B/lbm-°R 

5 

CP2(I) 

Specific heat - species 2 

B/lbm-°R 

5 

R1(I) 

Gas constant 1 - species 1 

B/lbm-°R 

5 

E2(I) 

Gas constant 2 - species 2 

B/lbm-°R 

5 

PVEMT(I) 

Relief valve vent pressure 

Psia 

5 

P(I) 

Initial tank pressure 

Psxa 

5 

TTMK(I) 

Initxal tank temperature 

■°R 

5 

TIN (I) 

Inlet gas temperature- species 2 

°R 

5 

TIME 

Time 

sec 

1 

VTANK(I) 

Tank volume 


5 

WOL(I) 

Initial mass of fluid 

Ibra 

5 

WD0TG(II,JJ) 

Inlet gas flowrate line II, node JJ 

Ibm/sec 

30,20 

WD0TG(II,1) 

Outlet gas flowrate line II, node JJ 

Ibm/sec 

30,20 

DQIN(I) 

Heat rate into tank 

Btu/sec 

5 
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5.1.12 TURBOE 

Turbonumn (Reference 1, Section 3.4.13) 



Symbol 

TURBOP(I) 

Description 

System component name for the turbopump I 

Dimensions 

5 

AREA! (I) 

Area of turbine inlet throat (Same as 
throat) turbopump I 

GG in^ 

6 

CEF1,CEF2(I) 

, Coefficients on turbine efficiency curve 


RPMD(I) 

Pump RPM at design speed 

rev/ min 5 

IPUMI(I) 

Line connected to pump inlet 


5 

IPUM0(I) 

Line connected to pump outlet 


5 

DIAD(I) 

Diameter pump design 

in 

5 

DIAP (I) 

Diameter pump actual 

in 

5 

DIAT (I) 

Diameter turbine 

in 

5 

gr(i) . 

Pump to turbine gear ration 


5 

ATR(I) 

Area turbine inlet 

. 2 
in 

5 

ATD (I) 

Area turbine discharge 

. 2 

in 

5 

KTD CD 

Discharge coefficient turbine 


5 

KTICD 

Discharge coefficient turbine inlet 


5 

3.1.13 PREG Pressure Regulator (Reference 1, Section 3. 

4.3.2) 


Symbol 

Description 

Units 

Dimensions 

PREG (I) 

System component name 


10 

PRIN(I) 

Line connected to regulator inlet 


10 

PROU (I) 

Line connected to regulator outlet 


10 

PRE (I) 

Regulator pressure discharge 

psia 

10 

3.1.14 Feedback Constraints (Reference 1. Section 3.4.18) 


Symbol 

Description 

Units 

Dimensions 

FEEDS (I) 

System component name 


6 

FBD(I) 

Line connected to feedback subroutine 
downstream (flow going into feedback) 


6 
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Symbol 

WFB(I) 

FBWC(I) 

FBPC(I) 

FBTC(I) 


Description 

Line connected upstream 
Flowrate feedback constraint equation 
Pressure feedback constraint equation 
Temperature feedback constraint equation 


Units . 


Dim,. 


0-T7S : 


6 

■f'6 

*6 


3.2 Output Description 

The output description is given in the output processor subroutine, 
Section 4.1. 

3.3. Typical Deck Set-up 

A typical ACPS model (Figure 1) was formulated to exercise all the 
program options. The hydrogen (H 2 ) side of the ACPS was modeled in detail, 
and the oxygen CO 2 ) side was simulated as constant pressure boundaries. The 
system consists of a liquid hydrogen (LH 2 ) storage tank feeding a turbopump 
assembly. The turbopurap increases the LH 2 pressure to 1421 psia and transmits 
the LH 2 to the heat exchanger which vaporizes the propellant. The heat 
exchanger and turbopiraip are driven by separate gas generators (gas-gas) which 
are operated in the feedback mode obtaining propellants from downstream of 
the heat exchanger. The heat exchanger also feeds an accumulator and three 
thrusters. Tlie following listing represents a typical set of card input 
and represents the schematic as shown in Figure 1, 
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N OMENCUT URE 

( I THRUSTOR 


COMPONENT 


PRESSURE 

boundary 


20* 1.7’' 
1 33 

SEC 


PI PL (JO) 
5‘ 1.7“ 


PRESSURE REGULATOR 


GAS GENERATOR 


GAS GENERATOR 

CHAM0 

2039"^R 

4.96 s£c 

450 PS) 


PIPL 

.. OA IBM 
5* 0.9576“ 


.PIPL@ 




PRESSURE BOUNDARY CONSTRAINT 


LIQUID LINE 
PIPL @ 

\ 1421 PS1 


’®VN 


PIPL ^ 

10' 6“ 

1347 PSIA JUNC10 


»2 PUMP 


TURBOPUMP 1 

V I V 


HEAT 

EXCHANGER 

HEATX0 


VALVG0 


)rURBOP0 


PIPL 

500°R 
1420 PSIA 
0.5' 6" 


^ PIPL (|^ 
5' 1 836“ 
PREG0 


P!PL@ \ 
10M.836" \ 
jSOOPSl \ 

T ^ PRE<K2) 

5' 1.836" 

7 PN 

r — ) 2 95" 

^ K PREG© 


LIQUID LINE 
PIPL © 


i/ 

HOT GAS 
^TURBINE V 


0.380 


GAS GENERATOR'' 

CHAM0 

1977®R 


y. 4 LDJVl 

0 21t5£^ FEEDS© p|pi^0 PIPL© 

-I PIPL© /'"’T 5' 400 PSl 550 PSl 

\ ; r 839" 20' 2 95" ' 

In 1 .^n LBM \„.«j LBM 


PRESSURE 

■^BOUNDARY 


5' 0 957" 


©XX® 

PIPL(4^ PIPL® 

0 r I 7"1 ^ 7" 

rTHRUSTOR(3^ 


/MR»4. \ 
500 PSl 
3000 LBP 


TANKOX 

® \ 
ACCUMULATOR} 
1346 P5J j 
. 50O°R / 

V H, y 


junoj(D 


20' I 7" 
, LBM 


PIPL(j2) 
ZO’ i 7” 
5.33 


VALVG© VALVG© 

i i 

PIPL© PIPL© 

0 I' 1 rj 0 r 1.7" 

|Thrustop01 

y^iR^4 \ 

/ 500 FSl \ 

pipu©X^200ufAp,p,@ 
20' 1.7" 20’ 1 7" 


VALVG© 


^VALVG© 


PIPL© PIPL© 

0 . 1 * 1 7" I |0 1 ' 1.7“ 

I THRU^OR®! 


15 99 LBM 
SEC 


' 510 PSl ^ 
3100 LBF 


Figure 1 


ACPS Demonstration Case 






SYSCOM 


5 'Al'lKC.2jPlPLLil6»PBL,3iPtFLi?7,PfPL»?6.TURBOP,l, 

PIPLLi17,pipi,,22»PBU»2«P5PL.23;CWAHi2»HEATX,1, 
PlPL,l8,'/ALV6«7|FlOi,,iq, ju iCLi 1 f P I PL ♦ 20 , 

PRFG,2»P!PL.2i,FFEDB,2,PfPL:,28iTANKD,lf 

PlPLi24|PREGi3,P!Fl<25,FEEQB*i; 

F I Pl » 1 1 PRtG 1 1 1 P I PL 1 2 i vJ ' JNC L 1 2 « P l PL » 3 , VA L VG > 1 8 P I PL « 4 , 

P!PL»5,VALVG,?,PIPL,6, 
„ . n « « , PIPL»7,VALVG,38PIPLi8f 

PbL 1 1 » P I PL * 9 8 JlJ.vcL » 3 1 R I FL 1 10 1 Val'-'G 8 4 1 P I PL * 11 8 CHA^ / 3 « 

P!FL«12eVAL'/Gf 5jPIPL» 138CHAM8 4f 
PIFL,14,VAL'/C',6-,PIPL,i5,CHAM,5, 


STaTf + Ft-’NcTlON'AL CCNSTRaIM3 

FUl\CT s C0^rT,i8C0^cT/2^P0'«C»i^ 

CGTAK’Kj 1 , 

CC^’CT,A,CO^CT,B^ 

CO^•CTs?8CO^CT,R,CO^CT^ 

CC'\cT 8 I J.cCKcT8ll8cONc’f’il2 8cO.XcT’»l^‘COMcTfl4,cnNcT8l5, 
TROCiI 1 

FBlCiliFeTC,?# FbPC81«FbPC<2 « 

STaTl 5 aRI'^'vJ'1«1»Pi031(aRJKj,182j Ctl58l6, PDELPf 1» 1« 1050 » » 

PG^2<^« l*399.*),p,i,i, 9?0i0CfPG«228 1 ;3 qo. , 

ARiKj,2ai08031»ARlNj',2i?,o;i?8 , 

ARINJai2iQi54j APEA7 ,381 j5,4 , 

ARn'’J»48l,Cij63, ARfNvj',4j2f 0;94aPEAT<4tl',5 4 , 
ARILji5<1»0«638ARINj', 5,2j 0 i 94 , ARgATj 5,l',5l4 , 

LENGTH, j,i,2000e 8 

TTFPP*22ii»500, «TTE>iP( 26 ,1,5 90 « 

sirultameous EQLATiOn, solvep 


K 3 20 « 

N s 20 « 

RAX IT = 80 ft 
»» 

7 Valve inputs 

*«• 

XK s 3R0, 9399, 3RC. 945, 2,0 
VAl - 3»B,7 8lO'si2»l4,ia ft 

VAo = 4i6, 8111,13/15,19 ft 

ft« 


ft* 


ftft three Cqi'BUsToFS + TWO CG«S 

ftft 


NCHAK 

ICHAF 

CV 

PCI 

TPRT(3} 
areat s 
NCQef 


s 5 ft 

<-1,1), 6, (5.a),l'3 % 

=/6/l2Rl,ft 


* 350 • >45o » , 2R5gn , ,B1 0 • * 
S 2R3C00, , 3l00, 

Cf228;2.7« 
a 4 ft 
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Junctions 

NLiNjU B 5|4,4 « 

IJUNC =/lO/(lil> il9, (i,2>#«20'. 

■'(3»i)»$t(3,2);wiQf{363)iel2{('5i4){«14ff 
WnOTG =/3C-^( 16» 1) •Bij(20»l)f2i70#^24^1)» e 211 » ( 2 ^ . i 5 1 1 » 00l5« 

if* lines 

XLENGL 15, ,15, .20. ,0.3 ,20e # 

0*i»2G,»C»l,2,0,5.0{ 

0.1,5,0,C.1,5.0,0.1» 

2Q,,,' oO,i5i;,15i,i5«( 

5,G,10. .5,C,15,,iQ, , 

5, 0,5,0, is, ■* 

E:^"IS c 28R0,07« 

EKlS(l6is 2R.0007« 

LIALO 5 6»l,3t05il,Rjlc6,lt8{ 
l,8ilt8,i»6,3,C5ili8, 

lQ?,d t8,1.8,i,8»l#S05 
6,l,6tl,^»l ,6,1,6,10, 

1,936,1.936,1,0576,6,1,1,936, 

1.936,1, £576,6, i» 

TM « 28R0,:5 
NCDeiLs 15R3, 2R2,11P3« 

NPIPL = 28 » 

GEE B 0,C« 
hNPR s 0.4* 

HNRE a o.e* 

hOEX s 0,25* 

TFIW s 6C0.’ 

NPLIN'F = 8P2i7ri,7P2,1>3R2i1,2^ 

^>DCTG =/3C/*'1i 15 »3»^975i I3ii) ,1,332-, ^5ii) il."5325i (7 il){i«3325, 
{9,11.15,^99, (10,15,5.33', (l?,t), 5. 33, <14 1 15,5,33 » 
TTeKiP =/3c/{ 9»i 5 , 5oO « I ( i 6, 1) 45 . , (27 f 1 5 , 5CQ. • { 23, 1 5 , 5 oO • ^ 

EGAS a 48,2813, 772, ’5* 

KA s 2 R 1.395 « 

«» 

«« PRESsupr pouneapies 

PGPBL * 550. ,5C0. ,4CQ, * 

NPBL a 3 « 

IPB a 9,23,27 * 
lEKD a 1,1,1 * 

«« 

Tank / accuhuaTcrs 

VTAnjK s 2C,,200, 

WOL = O.jlQCO, 

TUG a C. , 46, * 

CPl a .25,3,45^ 

CP 2 « ,25,3,45« 

R1 0 2R.993* 
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R2 a 2R.9?3« 

PVENT s 150Gi,60.a 
AVENT = 2P2,« 

DG IN z 2^0 I * 

LDMI s 2 fl,C« 

UDWO = C.i6« 

P(2) a 52. fl 

TTANK s 5aC.»45,« 

t>« HEAT E)iCH*NRER 

TyPF- = 2 » 

LI = 3 * 

LC = 1 * 

XI = 1 ,« 

XC = j 
TheC 550C. « 

HEOI =«»i7 » 

HEOO 5 16 tt 
HGII = 2 « 

RI = 3, « 

RC = 4.34 

Thick » o.C5 « 

* 

CEF2 » ~9,6« 

CEpi 5 5.2» 

RPMD = 500b0,» 

IPUHI s 16 ■ •» 

IPU^iO = 17 « 

GP ® 1 1 * 

DiAO =1, {} 

L I AP ” 1 1 * 

CUT = 24. Q o 
ATI * C.26 » 

ATD =l9» 4 

KTD =1, * 

KTI = 1. ft 
WTHCO^! = 3QP-'9,4t « 

«« PPPSSuRC f^-GijUTOR 
«« 

PR£ = 55C. i-OD, (4CDt* 

PRiN = 1»20«24 ft 
PROU - 2;2l«?-5 ft 

ftfl 

FEEPbaCK loops (PPESSUREf FLOV.R 4 TB F I XE{j TW£lEFORE Fb^C ^FbPC '^QT 

** PeOUJREd fBTC ONLY) 

FPI = 2T5I21 « 

FBD = 26f22 ft 

WDOTG = /'3o/{22,l)i2,70Qi{23>l5'.2t260, 0 . 169 , ( 26 , i ) , 0 i 2li «■ 

RFLAG = 7, ft 

ENDCAS 
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FSEOEDING PAGE BLANK NOT FILMED 


4.0 SUBROUTINES 

This section contains a list of the subroutines utilized in the SSACPS 
program and a detailed discussion on the 'more complicated subroutines. 
Descriptions of the remaining subroutines may be obtained from the program 
listing. 

Subroutine names and purpose follow; 


ADIAB 

BEGS 

BETA 

BDCINT 

BPROPG 

•BPEDPL 

CDDATA 

CHAM/TURB0P 

CHOICE 

CMON 

CWALL 

DP2R 

FCOMPl 

FC0MP2 

'FCOMP3 

FEEDB 

GATHER 

HEATEX 


Adiabatxc gas line 
Nonlinear equation solver 
Interpo lation 

Hollerith to integer conversion subroutine 

Calculates the gas properties at saturated conditions 

Calculates the liquid properties at saturated conditions 

Data Record Processing 

Bi-propellant combustor or turbopump 

Chooses heat exchanger configuration and calculation 
method 

Monopropellant combustor 
Thermal conductivity of wall 
Double precision to Real Conversion Function 
Heat transfer with friction 

Mach tiumber differential equation solved by Newton- 
Raphson iteration 

Isothermal pressure boundary differential equation 
by Newton-Raphson iteration 

Adiabatic Mach number differential equation solved 
by Newton-Raphson iteration 

Feedback constraints 

Gathers routines in order to create system as input 
by user. 

Heat exchanger calculations 
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HEATX 

Connects input parameters and line connections 
heat exchangers 

INITAL 

Tank initialization 

IN TEG 

Integration 

INTER2 

Interpolation 

INTEBP 

Interpolation 

ISOTH 

Isothermal gas line 

JUNCL 

Junction 

LEGS 

Single or double precision solution of an over- 
determined system of linear equations 

NUSLET 

Wusselt number 

OUTPRC 

Output processor 

PACK 

Character packing subroutine 

PEL 

Pressure boundary 

PEKFOR 

Computes dimensions 

PIPL/PIPLL 

Gas or liquid line with heat transfer and fret 

PREG 

Pressure regulator 

R2DP 

Real to double precision conversion junction 

SQZB 

Blank Dele tx on Subroutine 

SSPIPE 

Main program driver 

TAND 

Tank driver 

TANK 

T ank/ a c cumul at o r 

TRWLOD 

Numerical Data Input Subroutine 

VALVG 

Gas valve 

VALVL 

Liquid valve 

\vTRITE 

Writes output from heat exchanger 
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Symbo3 

BLK2 


PROPERTY SUBROUTINES 


BP BOP G 
BPROPL 

HPTCV 

HPTD 

HPTCP 

HPTTC 

HPTV 

H20H2 

LATENT 

PRDPTY 

P THE AT 

SAT 

THERM 

TPOGBl 

TPOCB2 


Description 

Initializes tables for use by PTHEAT 

Calculates saturated gas properties (oxygen and hydrogen) 

Calculates saturated liqxiid properties (oxygen and 
hydrogen) 

Computes the specific heat of hydrogen at constant 
volume 

Computes density of hydrogen from NBS data 

Computes specific heat for hydrogen at constant 
pressure 

Computes thermal conductivity for hydrogen from NBS 
data 

Computes viscosity of hydrogen from NBS data 

Calculates fluid properties of gas generator products 
(O^/H^ combustion) 

Calculates saturated gas and liquid enthalpies 
(hydrogen and oxygen) 

Thermodynamic property lookups (density, viscosity, 
heat capacity, thermal conductivity, Prandtl Number) 

Computes specific heat of hydrogen at constant 
pressure or volume (NBS data) 

Calculates the saturation temperature for hydrogen 
and oxygen 

Computes thermodynamic properties of oxygen 


Initializes Weber’s thermodynamic properties of oxygen 
Initializes Roder’s thermodynamic properties of oxygen 


Flowcharts 1 and 2 represent generalized flow charts of the SSACPS 
program logic showing the subroutine calling sequence. 
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Program Driver (SSPIPE) 



flowchart ACPS General Flowchart 
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Execute Subroutine calls 
as specif xed in the 
SYSCOM array 


Call system component 
models 


^re the\ 
'SYSCOM names 
\ legalV^ 


liIRITE 
ERROR 
Message 
No. 2 


STOP 

EXECUTION 


Set up functional values into 
array E as specified in the 
FUNCT array 


'^TJNCT variably 
\snames legaU-^'^ 


mTE 

ERROR 

Message 

No. 3 


WRITE 
FUNCTIONAL 
VALUES , 

f(i),J-=i.k/ 


EXIT 

GATHER 


Flowchart for Subroutine GATHER 
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Flov7chart 2 • 














Flowchart 3 . Flowchart for Subroutine OUXPFC 
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4. 1 OUTPRC 


Output Processor Subroutine 
Purpose 

The OUTPRC subroutine generates printed output at case termination for 
each of the xjser specified system component subroutines. 

Usage 

Galling Sequence 
0UTPRC 

Name 


Common Inputs 

Definition 

Units Dimensions 

AREAI(I,J) 

Inside line cross section area 
(line I, node J) 

. 2 
xn 

6,2 

ARE AT (I) 

Area of the combustor throat 

. 2 
in 

6 

ARINJ(I,J) 

Area of combustor injectors 
(line I, node J) 

2 

in 

6,2 

ATD(I) 

Area, turbine discharge 

. 2 

in 

5 

ATI (I) 

Area, Turbine Inlet 

. 2 

in 

5 

CDC(I) 

Discharge coefficient of the mono- 
propellant combustor injector 

— 

6 

CF(I) 
CP (I) 

Thrust coefficient 

Specific heat at constant pressure 

BTU 

6 

6 

lb-°R 

CSTAR(I) 

Characteristic exhaust velocity 

f t/sec 

6 

CV(I) 

Discharge coefficient of the bi- 
propellant combustor injector 


6 

CVEL(I) 

Isentropic spouting velocity 

ft/sec 

5 

DIALI (I) 

Inside diameter of line 

in 

30 

DIAT(I) 

Diameter of turbine 

in 

5 

DMVENT(I) 

Gas Flowrate through the vent 

Ibm/sec 

5 

UQIN(I) 

Heat into the accumulator/ tank 

BTU/sec 

5 

ETAT(I) 

Turbine efficiency 


5 
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OUTPRC 


Name 




Common Inputs 

Definition 

Units Dimensions 

FBL (I) 

Friction 

- 

30 

GR(I) 

Pump to turbine gear ratio 

- 

5 



BTU 

5 

HI (I) 

Heat* transfer coefficient, inside 
tube of the heat exchanger 

2 

in sec“R 

HCl(I) 

Heat transfer coefficient, outside 

BTU 

5 

2 


tube of the heat exchanger 

in sec“R 


HEAD (I) 

Equivalent heat transfer coefficient 
due to radiation 

BTU 

5 

sec°R 

ICHAM(I) 

Number of the line connected to the nth 
inlet of combustor (I) 

- 

6 

TEND (I) 

End to which pressure boundary (I) 
is connected 


6 

IPB(I) 

Used to identify which line is 


6 


connected to pressure boundary I 



IPUMI(I) 

Line connected to pump inlet 

- 

5 

IPHM0 (I) 

Line connected to pump outlet 

- 

5 



Ib^-sec 


ISP (I) 

Specific impulse of chamber (I) 

I 

6 

Ibm 

KTHER(I) 

Fluid conductivity line II 

BTU 

30 

hr-ft^F 

KTI(I) 

Total loss coefficient (liquid line) 

- 

30 

LAIIDA(I) 

Nozzle divergence coefficient 

-- 

6 

LDMI(I) 

Line connected to gas inlet to tank 

- 

5 

LDM(^(I) 

Line connected to gas outlet from tank 

- 

5 

LDWI 

Number of lines connected to liquid 
inlet to tank 

- 

5 

LDW0(I) 

Number of lines connected to liquid 
outlet to tank 

- 

5 

MACH (II, JJ) 

^ Mach number pipe 11 node JJ 

- 

20,30 

MEX(I) 

Chamber Mach number at exit 

- 

6 

MR(I) 

Mixture ratio 

- 

6 

MWC(I) 

Molecular weight of combustor combustion 
products 

lb 

lb mole 

6 
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OUTPRC 


Name 


Cotranon Inputs 

Definition 

Units 

Dimens ions 

NLINJU(I) 

Number of lines at junction (I) 

- 

10 

NODEL(I) 

Nodes in line I 

- 

30 

NPR(I,J) 

Prandtl number, line I, node J 

- 

30,2 

NRE(I,J) 

I^ynolds number of line I, node J 

- 

30,2 

NSYS 

Number system components 


1 

P(I) 

Pressure in tank I 

psia 

5 

PCI (I) 

Combustor chamber pressure (I) 

psia 

5 

PCN(I) 

Monopropellant combustor chamber 
pressure 

psia 

5 

PDELP(I) 

Pressure rise across turbopump 

psia 

5 

PDR0P(I) 

Pressure drop across injector 

psia 

5 

PE (I) 

Pressure at the exit of the 
combustor 

psia 

5 

PG(I,J) 

Pressure line I, node J 

psia 

30,20 

PGPBL(I) 

Pressure at pressure boundary 

psia 

6 

PGT(I,J) 

Total pressure line I, node J 

psia 

ft-lb^ 

30,20 

P0WP(I) 

Turbopump power 

r 

sec 

ft-lb 

6 

P0OT(I) 

Turbine power 

sec 

£ 6 

PPI(I) 

Pressure into turbopump 

psia 

5 

PP0(I) 

Pressure out of turbopump 

psia 

5 

PRIN(I) 

Line connected to pump inlet 

- 

5 

PROU(I) 

Line connected to pump outlet 

- 

5 

PTI (I) 

Turbine inlet pressure 

psia 

5 

Pl^jlCD 

Turbine outlet pressure 

psia 

5 

PW0(I) 

Coefficients used to describe the 
pump pressure rise/flowrate 
relationship for turbopump I 


4 
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Name 

Common Inputs 
R(I) 

BH0G(II,JJ) 

BH0L(I) 

EPMT(I) 

SYSC0M(I,J) 

TC(I) 

THRT(I) 

TTANK(I) 

TTEMP(II, JJ) 
TTI(I) 

TT0(I) 

TTTEMP(II,JJ) 

TWALL(II,JJ) 

U(II,JJ) 

UAD(II, JJ) 

VAD(I) 

VAI(I) 

7EL(II,JJ) 

VISC(I) 

MD0TG(II,JJ) 

WI(I) 

WN0Z(I) 

W0(I) 


Definition 

Units Dimensions 

Turbine velocity ratio 

- 

5 

Density of gas line II, node JJ 

Ibm/f 

30,20 

Density of liquid 

Ibm/ft^ 

30 

Turbine speed 

rev/min 

5 

System component names I component member - 

J comp onen t memb e r 

4,150 

Temperature of combustor 

“R 

6 

Combustor thrust 

Ibf 

6 

Initial temperature of accumulator/ 
tank 

“R 

5 

Fluid temperature line II, node JJ 

'’R 

30,20 ' 

Turbine temperature in 

®R 

5 

Turbine temperature out 

°R 

5 

Total temperature of fluid line II, 
node JJ 

°R 

30,20 

Temperature of the wall line II, node 

JJ °R 

30,2 

Overall heat transfer coefficient of 
liquid 

BTD 

hr-“R 

30,2 

Overall heat transfer coefficient of 
gas line II, node JJ 

BTU 

hr-^R 

30,2 

Line connected to valve discharge 

- 

20 

Line connected to valve inlet 

- 

20 

Velocity fluid line II, node JJ 

ft/sec 

30,20 

Viscosity fluid 

Ibm 

ft/sec 

30 

Flowrate fluid line II, node JJ 

Ib m/sec 

30,20 

Flowrate into junction 

Ibra/sec 

6 . 

Flowrate combustor nozzle 

Ibm/sec 

6 

Flowrate out of junction 

Ibm/ sec 

6 
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Name 

Common Inputs De finiti on Units Dimensions 


W0L(I) 

Weight of fluid in accumulator 

Ibm 

5 


Flowrate through piimp 

Ibra/sec 

5 

WT(I) 

Flowrate through turbine 

Ibm/sec 

5 

mHC0N(II,J) 

Fluid thermal conductivity, line II, 
node JJ 

BTO 

30,2 

hr-f t-®R 

XK(I) 

Discharge coefficient times the area 
of the valve 

. 2 

xn 

20 

XLENGL(I) 

Length of each section of pipe 

in 

1 


Subroutines Used 
WRITE 

Subroutine OUTPRC is shown in Flowchart 3, page -32 , 
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4.2 BDCINT 


Holier! th~tO“Integer Conversion Subroutine 
Purpose 

BCDINT computes an integer from an input array, in which each word 
consists of one left-justified hollerith character, 0 through 9 or a blanlc, 
with the remaining 5 characters being hollerith blanks. 


Usage 


Calling Sequence 

BCDINT NCOL, INT, N 

Input 

Argument Dimension 

NCOL N 


N 1 

Output 

INT 1 

N 1 

Restrictions 


Definition 

Array in which all words contain 
one left- justified BCD character. 
The remaining five characters of 
each word are blank 

The number of words in NCOL 

An integer corresponding to the 
BCD characters of NCOL 

The number of non-blank words in 
NCOL 


Each word of NCOL must contain one of the BCD characters 0 to 9 or a 
blank left- justified with the remaining 5 characters also blank. 

Subroutines Used 

SQZB 


Method 

BCDINT first calls SQZB to adjust the non-blank words of NCOL to the 
front. Then each word of NCOL is compared to each of the integers 0 through 
9 in BCD code, ^flien a match is found, the value to be returned (INT) is 
increased appropriately. If a non-valid character is found, the value of 
N is returned negative, in which case the magnitude of N indicates the 
digit in error. 

Subroutine BDCINT is shown in Flowchart 4, page 39. 
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Flowchart 4. Flowchart for Subroutine BCDINT 
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4,3 CDDATA 

Bata Be cord Processing Subroutine 
Purpose 

CDDATA reads an eighty character record in parts. The first 72 
characters may contain variable names, subscripts, data and other informa- 
tion xc^hile the remaining eight characters are available for sequencing or 
other user identification information* 

Usage 

CDDATA is called by TRWL(jSD for each input record. The record is 
read and resolved into the various fields, variable names, subscripts, 
control information, and individual data items. These fields are examined 
and categorized and loaded into a buffer (KBBP) with identifiers for use 
by TRT^0D. Numerous error checks are performed as the record is 
processed. 

Calling Sequence 
CDDATA , 

Input /Output 

All input and output is through common. 

Record Format 

The records acceptable to CDDATA are free form in that the 
information may be in columns 1-72. Columns 73-80 may be used 
for sequencing or for miscellaueoios information which is only read 
in and printed out on the tiser specified output unit. 

All columns betoeen two separators comprise a field. A 
separator may be an equal sign, a comma, an asterisk, a slash, 
a parenthesis, an R preceded by an integer, or apostrophes 
enclosing BCD characters. Blanks are permissible within all 
fields. An asterisk (*) except after an apostrophe for enclosing 
BCD characters terminate's the scan, and speeds processing of the 
record. 

Data information may be continued from one record to the next 
record by omitting the symbol name with storage of the information 
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continuing in the next consecutive cell. Data may be input 
for any location of a dimensioned variable (one-, too-, or 
three-dimensional array) by designating the subscript in 
parenthesis for a single dimensioned array (A = (3), 10.0), 
by designating the number of rows between slashes followed 
by the subscript for a too- dimensional array (B = /3/(3,3), 
10.9), or by designating the number of rows and columns 
between slashes followed by the subscript for a three- 
dimensional array (C = /3,2/ (2,2,3), 9.99). A three- 
dimensional array is the largest accepted by subroutine 
CDDATA. 

A BCD word preceding an equal sign (=) is a symbol, and 
all succeeding information is considered data for this symbol 
until the next symbol is encountered. Octal integers are 
input as 0123456700000, where the number must be preceded 
by the letter 0 and is stored from left to right as read. 

The values of all symbols are stored in the buffer in 
the same format in which they appear in the record, i.e. , 

BCD, integer, real, octal, or double precision. More than 
one symbol and its associated value (s) may appear in each 
record, i.e. X = 3.0, B = (2), 3.2, I = 4*. 

Error/Action Messages 

"ILLEGAL USE OF CEAMCTER IN OR AFTER COLUMN 

This message is written when an illegal character is located 
in any of the various fields. A few of the possible errors 
are as follows : 

a digit larger than 7 in an octal number, two 
periods in a real number, and a period in an 
exponent. 

Restrictions 

CDDATA will not store information in arrays having more 
than three dimensions. 

Subroutines /Functions Used 
BCDINT, PACK, SQZB 
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Method 

CDDATA reads a record using a 72A1, A6, A2 format, prints this record 
image if requested, and processes the various fields while storing the informa- 
tion in a buffer to be returned to Subroutine TRhT.OD. CDDATA scans the 
record until a field separator is located according to the separator of 
the field contents. The field and group type flags are defined as follows: 


Field Type 


Definition 


1 


Symbol follav7ed by an equal sign 
(6 characters or less) 


2 

3 

4 

5 

6 
7 


Single precision word (real) 

Single precision word (integer) 

Double precision word 
BCD string preceded by an apostrophe 
BCD string not preceded by an apqstrophe 
Octal integer 


Group Type - ’ Definition 

20 Dimension group 

30 Subscript group 

40 ‘ Repeat group 

The group type will always contain a field type 3, i.e, , 23, 33, or 43, 
representing a dimension, subscript, or repeat integer, respectively. 


T-lIhen the type flag has been determined and the field processed, the 
type flag follcwfed by the processed field is stored in the buffer. The 
procedure is repeated for all fields in the record until reaching an 
asterisk (*) or column 73, at which time a zero is stored in the next 
buffer cell signifying the end of the data. The number of buffer cells 
used for this record is stored in the buffer cell 109. The buffer has a 
storage capacity of 108 cells. 

Each symbol requires three buffer cells, one for the type flag, one 
for the field (symbol name), and one for the symbol storage address. Each 
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double precision number used three buffer cells* one for the type flag and 
tX'JO for the field (value). Each BCD field uses two buffer cells plus the 
number of buffer cells (six record columns per cell) needed for the field 
(BCD string); one for the type flag, one indicating the number of BCD words 
in the field, and the rest for the field. All the other fields only require 
two buffer cells, one for the type flag and one for the field. 

\-3hen a symbol is located, the table KSYMB is searched for this symbol, 
and if found, the address corresponding to the symbol is extracted from 
the table and stored in the buffer. If this symbol does not exist in the 
table KSYMB, it is an illegal input signified by a zero in the storage 
address cell of the buffer for this symbol. 


(1,2) 29.4, (2,3) 100.9, 1.2, 7.6* 
Definition 


Examples 

Data record Y = /5/ 

Buffer 

KEBE(l) = 1 
= Y 


= 23 
= 5 
= 33 
= 1 
= 33 
= 2 
= 2 
= 9.4 
= 33 
= 2 
= 33 
= 3 
= 2 

= 100.9 
= 2 
= 1.2 
= 2 


Type (symbol) 

Field (symbol name) 
Address of symbol 
Type (dimension integer) 
Field (number of raws) 
Type (subscript integer) 
Field (row subscript) 
Type (subscript integer) 
Field (column subscript) 
Type (real) 

Field (value of Y(l,2)) 
Type (subscript integer) 
Field (row subscript) 
Type (subscript integer) 
Field (column subscript) 
Type (real) 

Field (value of Y(2,3)) 
Type (real) 

Field (value of Y(3,3)) 
Type (real) 
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Buffer Definition 

KRBF(21) =7.6 Field (value of Y(4,3)) 

KRBF(22) = 0 Type (end of record) 


KRBF(109) = 22 Number of buffer cells used for 

this record 


Data record X = 123. , THIS IS NOT A SYMBOL, 012, 1.D5* 


KKBF(l) 


KRBF(17) 


Buffer 
= 1 
= X 


= 2 

= 123. 

= 6 

= 4 

= THIS I 
= S NOT 
= A SYMB 
= OL 
= 7 

= 120000000000g 

= 4 


Definition 

Type (symbol) 

Field (symbol name) 

Address of symbol 

Type (real) 

Field (value of X(l)) 

Type (BCD string not preceded by 
an apostrophe) 

Field (number of BCD words) 

Field (value of X(2)) 

Field (value of X(3)) 

Field (value of X(4)) 

Field (value of X(5)) 

Type (octal) 

Field (value of X(6)) 

Type (double precision) 



Field 1 

[ (value of X(7) and X(7)+l) 
Field 5 

Type (end of record) 


KKBF(109) = 17 


Number of buffer cells used for 
this record 


A functional flovjchart of this subroutine is shown in Flowchart 5, page 45. 




READ A RECORD 
AND PRINT 
ITS IMAGE 


SCAN A DATA 
FIELD AND 
DETERMINE ITS 
TYPE 


^THERE ANY^ 
ERRORS IN 
THE 

\^LD^ 


CONVERT THE 
FIELD 


WRITE ERROI 
MESSAGE 


STORE THE TYPE 
FLAG AND CONVERTED 
FIELD IN THE BUFFER 


WRITE A -1 
IN THE NEXT 
BUFFER CELL 


IS\. 

THIS TH^ 
END OF THE 
\ RECORD/ 


WRITE A ZERO 
IN NEXT BUFFER 
CELL AND THE 
NUMBER OF CELLS 
USED IN BUFFER 
CELL 109 


EXIT 

CDDATA 


Flowchart 5. 


Flowchart for Subroutine CDDATA 
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4.4 BP2R 


Double Precision to Real Conversion Function 
Purpose 

DP2R converts a double precision number stored in a type real variable 
to a real number. 


Functional Description 

DP2R calls the intrinsic function SNGL to perform the conversion 
from double precision to a real number. 

Usage 


Calling Sequence 
DP2R X 


Input 

Argument 

X 

Output 

Argument 

DP2R 


Dimension 

1 

Dimension 

1 


Error/Action Messages 
None 


Definxtion 

Value to be converted to real number 

Definition 
Converted value 


Restrictions 

None 
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4.5 R2DP 


Real to Double Precision Convert! on Function 
Purpose 

R2DP converts a real number to double precision to be stored in a 
type real variable. 

Functional Description 

R2DP calls the intrinsic function DELE to perform the conversion 
from real to double precision. 

Usage 

Calling Sequence 
X 


Dimension Definition 

1 Value to be converted to double 

precision 


Dimension Definition 

1 Converted value 

Error /Action Messages 
None 


R2DP 

Input 

Argument 

X 

Output 

Argument 

R2DP 
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4,6 PACK 


Character Packing Subroutine 
Purpose 

This subroutine packs BCD characters (one character per word) into 
array of six BCD characters per word. 

Usage 


Callxng Sequence 


PACK 


Input 

Argument 

NCABD 


N 

Output 

Argument 

LWORD 


N 


NCAPD, DWORD, N 


Dimension 


Definition 


K An array in which all words have 

one left-jiastified BCD character. 
The remaining fxve characters per 
word are blanks, 

1 , The number of words in NCARD 


Dimension 

Fn-i-sI 

y where 


Definition 

An array with six BCD characters, 
per word obtained by packing the 
input array NCARD ' . 


[] designated 
"greatest integer" 


1 


The number of words in DWORD 


Error/Action Messages 
None 


Restrictions 


None 

Sub routines /Functions Used 


None 



PACK 


Method 

PACK moves and combines BCD characters using bit manipulation and 
logical functions. The technique used to perform this function in the 
routine allows the subroutine to remain compatible for use on both the 
UNI VAC 1108 and the IBM 7094 computers. The follox^ing example shows a 
typical use of PACK: let NCAKD be the eight word array (one left- justified 

BCD character per word) "EXAMPLES" (i.e., NCAED = Ebbbbb, Xbbbbb.etc. ) , on 
output, DWORD (1) = EXAMPL and LW0RD(2) = ESbbbb. N should be input equal to 
8, and it will be returned equal to 2. NCAED is not restricted to alphabetic 
characters but may contain any alphanumeric characters. 

A flow chart indicating the logic used in Subroutine PACK is shown in 
Flowchart 6, page 50 . 
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Flowchart 6. Flowchart for Subroutxne PACK 
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None 


input may (in certain places), contain 
I processing* SQZB is called for this 


Definition 

An array in vzhich each word must 
he either all blanks or contain at 
least one noiX'-blanlc. BCD character 

The number of words in the array 
NCOL 


Definition 

If all the words of the NCOL array 
initially contained blanks or if all 
the words of the NCOL array contained 
at least one non-blank BCD character, 
then the NCOL array is unchanged. If, 
however, there is a mxxture of the 
then the output NCOL array will have 
non-blank words at the beginning of 
the array and the remaining words of 
the array filled with blanks. See 
example in Method. 

The number of non-blank words in the 
array NCOL 
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Restrictions 

None 

Sub routines /Functions Used 
None 

Method 

SQZB removes blanks from unpacked (one character per word) arrays of 
variable length* Non-blank words are left-justified xn the array and the 
remaining words are blanked. 

Example 

Let the input array NGOL be ten words long with each word containing 
one BCD character (left- justified with blank fill): b, b, S, A, M, b, P, L, 
b, E (where ^'b‘* represents a blank)* On return from SQZB the array NCOL 
will be: S, A, M, P, L, E^ b, b, b, b, one character per word, with N=6* 

A functional flowchart of this subroutine is shown in Flowchart 7, page 53 . 
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Flowchart 7 ♦ 



Flowchart of Subroutine SQZB 
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4.8 TRWLOD 


Numerical Data Input Subroutxne 
Purpose 

TRifLOD is an input control routine which takes data contained 
by a buffer and stores the value (s) of each variable in specified core 
location(s) . 


Usage 

Calling Sequence 

TR^TLOD I IN, I0UT, K 


Input 

Argiiment 

IIN 

I0UT 


Dimension 

1 

1 


Definition 

Input BCD tape unit 

0utput unit for record images. 
Suppress output if equal to 0. 


Output 

Argument Dimension 

K 1 


Definition 

Error flag = 2, Indicates an 
error 


Error/Action Messages 

1) "ERROR IN SUBSCRIPT (SEE BUFFER WORD ^)" 

This error message is written for any of the following sub- 
script errors. The type flag from the buffer designates 
that a dimension integer is specified when the dimension 
has already been defined, the column subscript is being 
used without the number of rows having been specified, or 
the third subscript is being used without the number of 
columns being specified. 

2) "INCORRECT DATA TYPE. NO CONVERSION POSSIBLE. 

WAS EXPECTED, FOUND ON CARD (SEE BUFFER WORD . )" 

This message is written when the type of the given value 
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is not consistent with the specifications in bits 18-20 
of the address word NADR and the value cannot be converted 
without changing its value. 

3) " IS AN ILLEGAL SYMBOL" 

A 

this message is written when the variable name in the 
buffer is not a legal symbol. 

4) "ILLEGAL DATA CAED CONTAINED IN BUFFER" 

This message is written when the type flag is greater than 
7, a symbol has not been defined (J=0) and the first buffer 
cell contains a value other than a 1; i.e. , the first data 
card is a continuation card, or when non-integer subscript, 
dimension, or repetition values have been used. 

5) "ERROR INVOLVING REPETITION OF DATA FIELD (SEE BUFFER WORD )" 

This message is written when either a negative or zero 

number of repetitions have been requested or repetition of 
a non- repeatable field is indicated; i.e. , symbol name, 
subscript, etc. 

Subroutines Used 

CDDATA 

Functions Used 

DP2R, R2DP 

Method 

TR^-ILOD calls Subroutine CDDATA to process a data record and returns 
the information contained by a buffer. TRWLOD processes this information 
storing the value (s) in the location (s) specified according to the given 
subscript (s) (if any). TRWLOD checks the type of the value in the buffer 
with the type specified for each variable, and if the types are not com- 
patible, the value is converted providing the conversion can be accomplished 
without altering the values. For instance, integers will be converted to 
real numbers and real numbers converted to integers when conversion does 
not cause truncation of the fractional part. 
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If the type flag from the buffer designates that a symbol is in the 
next buffer location, the symbol is stored in the variable "NAME." The 
buffer cell (stored in "NADR”) following the symbol has the storage address 
in bits 21-35 and the type (real, integer, BCD, octal, etc.) of the vaiue(s) 
expected in bits 18-20 (stored in "LL"). If the storage address is zero 
and the symbol is not "ENDCAS," this is an illegal symbol, and error message 
number 3 is written. If the symbol was "ENDCAS," control is returned to 
the calling subroutine. 

A type flag equal to a too, three, four, five, six, or seven (type 
real, integer, double precision, BCD characters limited by apostrophes, 

BCD characters, or octal, respectively) designates the type of the value 
in the next buffer cell. If this type does not agree with "LL," the value 
is converted, if possible, without changing its value, e. g., integer to 
real; othert-rise, error message number 2 is written. For any value being 

converted, the message "INPUT VALUE CONVERTED TO is written. 

If the type flag is the same as the value of "LL," "LL" xs equal to zero, 
or if the value has been converted, the value is stored according to its 
subscripts (if any). A type flag equal to 1 causes error message number 4 
to be wrxtten. 

A type flag of 33 or 23 (type subscript integer or dimension integer, 
respectively) designates a subscript or dimension information is contained 
in the next buffer cell. T^hen receiving subscrxpt or dimension information, 
the subscript or dimension is determined, and if not correctly specified, 
error message number 1 is written. 

A type flag of 43 (repeat integer) signified repetition of the follow- 
ing field is to take place the number of times specified in the next buffer 
location (location following the 43). 

If "LL" is equal to zero, the values are stored as they are input 
xjithout any conversion. A type flag equal to zero in the buffer indicates 
that all the information in the buffer has been processed. 

After any error message is written, the buffer is dumped for use in 
diagnosing the error. 

A functional flowchart of this subroutine is sho™ in Flowchart 8, page 57 
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ENTER 

TMvfLOD 


WRITE ERROR 
MESSAGE 5 


^ IS \ 
THE TYPE ^ 
GREATER THAN 
s. 7 ^ 


WRITE ERROR 
MESSAGE 4 


Flowchart 8. Flowchart for Subroutine TRWLOD 











4.9 PIPL/PIPLL 

Fluid line subroutine simulating heat transfer and frictron. 

Purpose 

Subroutine PIPL calculates the gas line properties and PIPLL 

calculates liquid line fluid properties that result from flowing through a 

line having heat transfer and friction. 

Usage 

Calling sequence: 

PIPL 11 

Input : 


Argument 

Dimension 

Definition 


II 

1 

Gas line number 


Entry Point 




PIPL 

II 



Output: 




Argument 

Dimension 

Definition 


II 

1 

Liquid line number 


Name Common 
Variables 

Dimension 

Name Common 
Variables 

Dimension 

ALPHA 

30 

FRL 

30,20 

AREAI 

30 

GC 

1 

AREAK 

30 

GEE 

1 

AREA0 

30 

HI 

30,2 

CP 

30 

HHPR 

1 

DIALI 

30 

HHRE 

1 

DIAL0 

30 

H0 

30,2 

EMIS 

30 

• H0EX 

1 

FB 

1 

HR 

5 
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Name Common 
Variables 

Dimension 

Name Common 
Variables 

Dimension 

HRAD 

30,2 

PI 

1 

HRS 

5 

RFLAG 

1 

KA 

2 

RGAS 

2 

KLT 

5 

RH0G 

30,20 

KIHER 

30 

RH0L 

3 

LENGTH 

5 

TENV 

1 

MACH 

30,20 

TTEMP 

30,20 

NGR 

30,2 

TTTEMP 

30,20 

N{5DEL 

20 

TWALL 

30,2 

NPLINE 

30 

UA0 

30,2 

NPR 

30,2 

VEL 

30,2 

NRE 

30,2 

Vise 

30 

PG 

30,20 

TO0TG 

30,20 

PGT 

30,20 

mHC0N 

30 

Sub routines /Functions 

used 



0PTCP, 0PTCV, 0PTTC, 

0PTV, HPTCP, 

HPTCV, HPTTC, HPTV, 0PTD, 

HPTD, FC0MP1 


Flowchart 9, page 61, shows the logic used in Subroutines PiPL/PIPLL. 
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PERFORM ENERGY BALANCE ! 
AND CALL FC0MP1 FOR MACH 
NUMBER SOLUTION 


COMPUTE DO^aJNSTREAM 
FLUID PROPERTIES 


Flowchart 9 . 


Flowchart for Subroutine PIPL/PIPLL 
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4.10 HEATEX 


Heat Exchanger Subroutine 


Purpose 

This subroutine calculates the heat exchanger outlet conditions given 


the fluid inlet conditions. The heat exchanger calculations are automatically 

performed at twenty (20) segments along the heat exchanger. 

Usage 

Calling Sequence 

HEATEX LN,LI ,L0 ,TI1 ,T01 ,TI2 ,T02 ,PI1,P01,PI2 ,P02 ,WDI ,1^00 ,QT,XI1 ,X01 


Input 

Argument 

LN 

LI 

-L0 

Til 

T01 

PIl 

P01 

TOI 

Ti/D0 

XII 

X01 

0utput 

Argument 

TI2 

T02 

PI2 

P02 

QT 


Dimension 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 


Dimension 

1 

1 

1 

1 

1 


62 



HEATEX 


Name Common 
Variables 

Dimension 

Name Common 
Variables 

Dimension 

ARFL 

5 

PHIl 

1 

CPI 

21 

PHI 2 

1 

CP0 

21 

PHIS 

1 ' 

DHY0 

5 

PI 

21 

LIST 

21 

P0 

21 

DSI 

5 

PEI 

21 

DTB 

5 

PR0 

21 

FD 

5 

Q 

21 

FN 

5 

QUALI 

21 

HE 

21 

QUAL0 

21 

HI 

21 

REI 

21 

H0 

21 

RE0 

21 

HH 

21 

RI 

5 

KI 

21 

R0 

5 

K0 

21 

SUI 

21 

KW 

21 

SV0 

21 

LENGTH 

5 

THICK 

5 

LT 

5 

TI 

21 

MUI 

21 

T0 

21 

MU0 

21 

TWI 

21 

NT 

5 

TW0 

21 

P 

5 

TYPE 

5 


Subroutxnes/Functions Used 

PRl^PTY, SAT, LATENT, BPR0PL, BPR0PG, NUSLET, WALL 


Flowchart 10, page 64, shows the logic used in HEATEX. 
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Flowchart 10. Flowchart for Subroutine HEATEX 
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4.11 ADIAB 


Adiabatic Gas Line Subroutine 
Purpose 

TTie purpose of this subroutine is to calculate fluid properties in a 
Ixne where no heat is transferred. 

Usage 

Calling Sequence 
ADIAB II 

Input 

Argument Dimension 

II 1 


CP 

30 

DIALI 

30 

FB 

1 

FEL 

30, 20 

GC 

1 

KA 

2 

MACH 

30, 20 

Name Common 


Variables 

Dimension 

NPLINE 

30 

NEE 

30, 2 

PG 

30, 20 

PGT 

30, 20 

RFLAG 

1 

EGAS 

2 

EH0G 

30, 20 

TTEMP 

30, 20 

TTTEMP 

30, 20 

VEL 

30, 2 

Vise 

30 

WD0TG 

30, 20 

XLENGTH 

30 
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ADIAB 


Sub routines /Fmctlons Used 

^PTCP, 0PTCV, 0PTV, HPTCP, HPTCV, HPTV, PC0MP3 

Plowchart 11 j page 67, shows the logic used in Subroutine ADIAB. 


66 









4. 12 CIKDN 


Monopropellant Combustion Subroutine 
Purpose 

The purpose of this routine is to calculate the combustion parameter 
for a monopropellant conijustor assuming a constant ratio of specific heats. 

Usage 

Calling Sequence 
CMON I 

Input 

Argument Dimension 

I 1 


Name Common 
Variables 

Dimension 

Name Common 
Variables 

Dimension 

AT 

1 

N0DEL 

20 

AKEAE 

6 

NPLIHE 

30 

ARE AT 

6 

PATH 

1 

AEINJ 

6, 2 

PCI 

6 

CDC 

6 

PCN 

6, 2 

CTW 

1 

PDR0P 

6 

CV 

6, 2 

PE 

6 

GC 

1 

PG 

30 ' 20 

I CHAM 

6, 2 

PGT 

30, 20 

ICM0N 

1 

RGAS 

2 

ISP 

6 

RELAG 

1 

ISPI 

6 

THRT 

6 

KA 

2 

TTEMP 

30, 20 

LAMDA 

6 

TITEliP 

30, 20 

HEX 

6 

WD0TG 

30, 20 

NCHAM 

i 

WN0Z 

6 


68 



QION 


Subroutines/ Functions Used 
OPTCP, OPTCV, HPTCP, HPTCV 

Flowchart 12, page 70, shows the logic used in Subroutine CMON. 


69 





CALCULATE CONSTRAINT 
EQUATION FOR CHAMBER 
PRESSURE 

is 





CALCULATE CONSTRAINT 
EQUATION FOR SPECIFIC 
IMPULSE 



CALCULATE CONSTRAINT 
EQUATION FOR THRUST 


RFLAG 


WRITE ROUTINE 
PRINTOUT 


Flowchart 12 • 


Flowchart for Subroutine CM0N 










4.13 CHAM/TURB0P 


Bi-propellant Combustor and Turbopump Subroutines 
Purpose 

The bi-propellant combustor models either a LO 2 /H 2 gas generator or a 
LO 2 /H 2 combustor. The combiostor performance data is calculated based on 
curve fit data input as a function of mixture ratio for the propellants 
combining xnto the combustor. The combustor can operate as a gas generator 
and therefore be connected to either the turbopump or the heat exchanger 
model. The turbopump model simulates a pump connected through a gear box 
driven by a gas turbine. 

Usage 


Calling Sequence 




CHAM I 




Input 




Argument 

Dimension 



I 

1 



Entry Point 




TURB0P I 




Input 




Argument 

Dinsnsion 



I 

1 



Name Common 


Name Common 


Variables 

Dimension 

Variables 

Dimension 

AJ 

1 

C0NCT 

3,6 

AKEAT 

6 

CP 

30 

AEINJ 

6,2 

CPCX 

6 

ATD 

5 

CP2 

5 

ATI 

5 

CSTAR 

6 

CEFl 

5 

CS2 

4 

CEF2 

5 

CV 

6,2 

CF 

6 

CVEL 

1 


71 



CHAM/TURB0P 


Name Common 
Variables 

Dimension 

Name Common 
Variables 

Dimension 

DIAD 

- 5 

P0OT 

' 5 

DIAP 

5 

P0¥1 

5 

DIAT 

5 

P0N2 

5 

ETAT 

1 

P0W3 

5 

EB 

1 

- PPI 

5 

GC 

1 

PP0 

5 

GR 

5 

PTI 

5 

I CHAM 

6,2 

PW0 

5 

IPUMI 

5 

PW0 

5 

IPUM0 

5 

PWl 

5 

ISP 

6 

PW2 

5 

ISPT 

6 

R 

1 

IS2 

4 

BFLAG 

1 

KA 

2 

RGAS 

2 

KTD 

5 

EH0L 

3 

KTI 

5 

EPMD 

5 

MR 

6 

RPNT 

5 

MWC 

6 

TC ' 


Mt-72 

4 

TC2 

4 

NC0EF 

1 

THRT 

6 

N0DEL 

20 

TTEMP 

30, 20 

NPLINE 

30 

TTI 

5 

PCI 

6 

TT0 

5 

PDELP 

5 

TTTEMP 

30, 20 

PDR0P 

' 6 

U 

1 

PG 

30,20 

WD0TG 

30, 20 

PGT 

30,20 

WN0Z 

6 

PMR 

6 

WP 

5 

P0WC 

5 

WT 

5 

P0WO 

P0WP 

5 

5 

TOGO 

5 
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CHAiI/TURB0P 


Sub routines /Functions Us0d 
0FTCP, 0PTCV, HPTCP, HPTCV 

Flowchart 13, page 74, shox^s the logic used in Subroutine CHAM/TURB0P. 
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COMPUTE WN0Z, MR, ISPT, 
CSTAR, PMR, MtTC, TC, CPCX 
AND CF 



CALCULATE CONSTRAINT EQUATIONS 
FOR THE INJECTOR OXYGEN FLOW 
RATE AND THE INJECTOR HYDROGEN 
FLOW RATE 



Flowchart 13. Flowchart for Subroutine CHAM 
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A 


CALCULATE CONSTRAINT 
EQUATION FOR CEAKBER 

PRESSURE 


'RFLAG 


NGGTP 


EXIT 

CHAM 


COMPUTE PUMP AND TURBINE 
PARAMETERS PPI, TO, PP0, 
KPMP, P0OT, RPMT, TTI, PTI, 
CVEL, R, U,' ETAT, P0WT, TT0, 
WT, AND PT0 _ 


CALCULATE CONSTRAINT EQUATIONS 
FOR THE TURBINE GG GAS FLOW RATE 
AND THE PUMP-TURBINE POWER 


SETUP PIM> discharge 
PARAMETERS WD0TG, TTEMP, 
AND PG 


WRITE 

ROUTINE 

PRINTOUT 


RFLAG 


Flowchart 13. Flowchart for Subroutine CHAM (Continued) 




ISOTH 


4.14 IS0TH 

Isothermal Gas Line Subroutine 
Purpose 

This subroutine simulates an isothermal gas line and determines downstream 
fluid conditions. 

Usage 

Calling Sequence 
IS0TH II 

Input 

Argument Uimension 

II 1 


Name Common 
Variables 

' Dimension 

Name Common 
Variables 

Dimension 

CP 

30 

PGT 

30,20 

DIALI 

30 

RFLAG 

1 

FB 

1 

RGAS 

2 

FRL 

30,20 

RH0G 

30,20 

GC 

1 

TTEZ'IP 

30,20 

KA 

2 

TTTEMP 

30,20 

MACH 

30,20 

VEL 

30,2 

NPLINE 

30 

Vise 

'30 

NRE 

30,2 

WD0TG 

30,20 

PG 

30,20 

XLENGL 

30 


Sub routines /Functions Used 


0PTCP, 0PTLP, 0PTV, HPTCP, HPTCV, HPTV, FC01®2 

Flowchart 14, page 77 , shows the logic used for the Subroutine ISOTH. 
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CALCULATE FLUID PARAMETERS 
RGAS, KA, VEL, TTTEMP, RH0G, 
NRE, AND FRL 


1 


CALL FC0MP2 FOR STATIC 
Pips SURE SOLUTION 



Flowchart 14. Flowchart for Subroutine IS0TH 
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PRECEDING PAGE BLANK NOT FILMID 


SYMBOLS 

ARINJ 

ALPHA 

AREAI 

AREAT(I) 

AREG(I) 

ATD(I) 

AVEKT(I) 

CDC(I,J) 

CEFI(I) 

CEF2(I) 

CP 

CPO, CPOG, CPI 

CPI, 2(1) 

CSTAR (I) 

CV 

CV(I,J) 

CVEL(I) 

DELXL(II) 

DHYO 


5.0 NOMENCLATURE 

DESCRIPTION 
Oxidizser area injector 

Angle line makes with horizontal 

Inside area of line 

Combustor throat area combustor I 

Maximum area of regulator asstiming 
sonic flow of regulator (I) 

Turbine duct diaxhage area 

Flow area vent valve tank (I) 

Discharge coefficient through injector 
monoprop ellant 

Turbine efficiency coefficient curve 

Turbine efficiency coefficient curve 

Specific heat at constant pressure 

Specific heat at constant pressure 
outside, gas, inside in heat exchanger 

Specific heat at constant pressure of tank 

Characteristic exhaust velocity combustor ! 

Specific heat at constant volume of tank I 

Discharge coefficient through injector, 
bipropellant 

Isentropic spouting velocity of turbopump 
Line segment (line II) 

Hydraulic diameter of heat exchanger 


DIAD(I) Pump design diameter of turbopump (I) 

DIALI, DIAL0(II) Diameter line, inside, outside (line II) 


UNITS 

2 

in 

deg 

. 2 

xn 


. 2 

xn 

, 2 

in 


BTU/ibm-°R 

BTU/lbm-°R 

I BTU/lbm-‘’R 
I ft/sec 

BTU/lbm-°R 

(I) £t/sec 
ft 

Ibm/ in- 
sec 

in 
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SYMBOLS 

DESCRIPTION 

UNITS 

DIAP(I) 

Puinp diameter of turbopump (I) 

in 

DIAT(I) 

Turbine rotor diameter 

in 

DX 

Length of heat exchanger 

in 

EMIS 

Eraisslvity 

- 

ENTH,0,I,G,L 

Enthalpy outside, inside, gas, liquid 
in heat exchanger 

BTU/lbm 

EP(I) 

Pump efficiency of turbopump (I) 


ETAT(I) 

Turbine efficiency 

- 

FBPC(I) 

Pressure constraint of junction. (I) 

- 

FBWG(I) 1 

Flowrate constraint of junction (I) 

- 

FFACTI, 0 

Friction factor inside, outside in 
heat exchanger 

- 

FRL 

Friction factor 

- 

GAM 

Specific heat ratio C /C, 

* p V 

- 

GC 

32.174 

Ibm-f t/ 
sec^lbf 

HE 

Equivalent heat transfer coefficient 
in heat exchanger 

BTD/in^- 
°R -hr 

HI 

Inside film coefficient 

BTU/hr- 

ft^-°R 

HI, 0 

Heat transfer coefficient inside, 
outside in heat exchanger 

BTU/in^- 

sec-°R 

H0 

Outside film coefficient 

BTU/ft^- 
°R -br 

HR 

Hydraulic Resistance 

sec^/f t- 



, z 

in 

HRAD 

Radiation heat transfer coefficient 

BTU/hr- 

ft^‘« 
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ISP (I) 

K 

KA 

KAY 

KLT 

KTD(I) 

KTHER 

KTI (I) 

MACH (II, JJ) 
MARTO(I) 

MR(I) 

MUI 

NCR 

NLIKJU(I) 

NPR 

NRE 

P 

P 

PB(I) 

PCI(l) 
PCORI, 0 

PDROP(I,J) 


Vacuum specific impulse combustor I 

Fluid thermal conductivity in heat 
exchanger 

Specific heat ratio 

Ratio of specific heats combustor 

Total loss coefficients; bends, contractions, 
etc. 

Turbine discharge coefficients of 
turbopump (I) 

Fluid thermal conductivity 

Turbine inlet coefficient of turbopump(I) 
Mach niffiiber down pipe, line II, node JJ 
Martinelli parameter 

Mixture ratio combustor I 
Viscosity fluid in heat exchanger 

Grashof number 

Number of inlet or outlet lines at 
junction (I) 

Prandtl number 

Reynolds number 

Fluid pressure in heat exchanger 
Pressure in tank I 

Pressure boundary constraint equation 
in the Ith pressure boundary 

Chamber pressure, combustor I 

Pressure drop correction for two phase 
flow, inside, outside in heat exchanger 

Pressure drop through injector J, 
combustor I 


sec 

BTU/ in- 
sec- R 


BTU/sec- 

ft-°R 


Ibm/sec 

Ibm/ in- 
sec 


psia 

psia 

psia 

psia 
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PG(II,JJ) 

Pressure, line II, node JJ 

psia 

PG(II,JJ) 

Pressure, line II, node JJ 

psia 

PHIl 

Correction for two phase flow boiling 
in heat exchanger 

- 

PHI2 

Entrance effects correction in heat 
exchanger 

- 

PHI3 

Asynunetric heating correction in heat 
exchanger 

- 

PHI4 

Surface roughness correction in heat 
exchanger 


P0W,0,1,2,3 (I) 

Coefficients in pump power design curve 

- 

POW(I) 

Pump power of turbopump (I) 

ft-lbf/ 

sec 

POWT(I) 

Turbine power of turbopurap (I) 

ft-lbf/ 

sec 

PRE(I) 

- Pressure regulator doVJnstream of 
regulator I 

psia 

PWO,l»2, (I) 

Coefficients in pump pressure design 
curve 


Q 

Heat flux in heat exchanger 

ETU/sec*“ 

. 2 

xn 

R(I) 

Gas constant in tank I 

B/lb-°R 

RGAS(I) 

Gas constant in line 

f t-lbm 
lbf-°R 

RHOG,RHOL(II,JJ) 

Fluid density, gas, liquid 
(line II, node JJ) 

Ibm/f t^ 

RHOL(I) 

Density in tank I 

Ibm/ft^ 

RI, RO 

Heat exchanger radius 

in 

RPMD(I) 

Pump design speed 

RPM 

EPMP(I) 

Pump speed of turbopump (I) 

EPM 

RPMT(I) 

Turbine Speed 

RPM 

SVO, 0,I,G,L 

Specific volume outside, inside, gas, 
liquid in heat exchanger 

in^/lbm 
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SVBAR 

Change in specific volume average in 
heat exchanger 

ix? / Ibm 

T 

Average gas temperature in tank 


TC(I) 

Chamber temperature combustor I 

°R 

TENV 

Temperature environment 

°R 

THICK 

Wall thickness xn heat exchanger 

in 

THRTd) 

Thrust, combustor I 

Ibf 

TIME 

Time 

Sec 

TIN (I) 

Gas inlet temperature in tank I 

°R 

TPI (I) 

Pump inlet temperature of turbopump (I) 

°R 

TP0(I) 

Pump discharge temperature of turbopump 

(I) 

°R 

TTEIS 

Turbine isen tropic discharge temperature 
of turbopump (I) 

°R 

TTEMP(II,JJ) 

Static temperature fluid line II, 
node JJ 

°R 

TTI (I) 

Turbine inlet temperature of turbopump (I) 

°R 

TTO (I) 

Turbine discharge temperature of turbo~ 
pump (I) 

°R 

TTTEMP(II,JJ) 

Total temperature stagnation fluid, line 
II, node JJ 

“R 

TWI,TWO 

Discharge temperature in heat exchanger, 
inside, outside 

°R 

Tt^ALL 

Wall temperature 

°R 

UA0 

Overall heat transfer coefficient 

BTU 

°R sec 

VEL(II,JJ) 

Velocity fluid, line II 
node JJ 

ft/sec 

Vise 

Viscosity 

Ibm/f t- 
sec 

VTANK(I) 

Volume of tank I 

ft^ 
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w 

W(I,2) 

T-7D, WD0TG(II,JJ) 
WGASI, WGAS0 
WN0Z(I) 

WOL(I) 

WP (I) 

WT (I) 

XK (1) 

a 

0.1714 X 10^ 


Flowrate Ibm/sec 

Oxidizer flowrate combustor I Ibm/sec 

Fuel flowrate combustor I Ibm/sec 

Flowrate, line II, node JJ Ibm/sec 

Floxnrate xn heat exchanger, inside, outside Ibm/sec 

Nozzle flowrate combustor I Ifam/sec 

Mass of liquid in tank (I) Ibm 

Pump flowrate Ibm/sec 

Turbxne flovTrate Ibm/sec 

2 

Cross sectional area of valve times in 

valve discharge coefficient of valve I 


Mass ratio, m^/ (m^ + m^) 
Stefan- Boltzmann constant 
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